A New Phase of Tethered Membranes: Tubules 
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We show that fluctuating tethered membranes with any intrinsic anisotropy unavoidably exhibit a 
new phase between the previously predicted "flat" and "crumpled" phases, in high spatial dimensions 
d where the crumpled phase exists. In this new "tubule" phase, the membrane is crumpled in one 
direction but extended nearly straight in the other. Its average thickness is Eg ~ i'^' with L the 
intrinsic size of the membrane. This phase is more likely to persist down to d = 3 than the crumpled 
phase. In Flory theory, the universal exponent i/t = 3/4, which we conjecture is an exact result. We 
study the elasticity and fluctuations of the tubule state, and the transitions into it. 
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Tethered mombranos are of great interest in large part 
because their behavior is much richer than that of poly- 
mers, their one-dimensional analog. Specifically, poly- 
merized membranes have been predicted [1] to undergo 
a "crumpling" transition between the "crumpled" and 
long-ranged orientationally ordered "flat" phases. This 
apparent violation of the Mermin- Wagner theorem is 
made possible by "anomalous elasticity" [1], [2]: ther- 
mal fluctuations inflnitely enhance the membrane's effec- 
tive bending rigidity k, stabilizing the orientational order 
against these very fluctuations. 

Most past theoretical work [3] has been restricted 
to isotropic membranes. Here we consider intrinsically 
anisotropic membranes and find that this seemingly in- 
nocuous modification has profound and surprising con- 
sequences: an entire new and heretofore unanticipated 
phase of the membrane, which we call the "tubule" phase 
ubiquitously intervenes between the crumpled and "fiat" 
phases (see Fig.l). Only in the special case of perfectly 
isotropic membranes, which follow a path like P2, is a 
direct crumpled-to-fiat (CF) transition possible. Generic 
paths like Pi have two phase transitions, crumpled-to- 
tubule (CT), and tubule-to-fiat (TF), which we are cur- 
rently studying in an e-expansion [4]. 

There are a number of possible experimental realiza- 
tions of anisotropic membranes. One is polymerized 
membranes with in-plane tilt order. Fluid membranes 
with such order have already been found [5]; it should 
be possible to polymerize these without destroying the 
tilt order. Secondly, membranes could be fabricated by 
crosslinking DNA molecules trapped in a fluid membrane 
[5]. Performing crosslinking in an applied electric fleld 
would align the DNA and "freeze in" the anisotropy in- 
duced by the electric field, which could then be removed. 

Simulations could be done on, e.g., triangular or rect- 
angular nets of balls and springs with all of the spring 



constants in one direction different, by a factor of order 
2 or so, from those in the other direction. Equivalently, 
one could have different bond lengths in the two direc- 
tions, or use second nearest neighbor springs of differ- 
ent strengths to create different bend stiffnesses in the 
two directions. Any such modification whatsoever will 
lead, upon renormalization, to a membrane with all of 
the anisotropic terms we consider here, and, hence, will 
fafl into the universality class of our model. 




FIG. 1. Phase diagram for anisotropic tethered membranes 
showing the flat, crumpled and new tubule phases. 

The defining property of the tubule phase is that the 
membrane is crumpled in one direction (y) , but "flat" 
in the other. Its average shape is a long, thin cylinder of 
length Ry = L°xO{l) and radius Rg{L±) -C L±, where 
Ly and L± are the intrinsic dimensions of the membrane. 
The tubule radius Rq, and its undulations hrms trans- 
verse to the y-axis, obey the scaHng laws: 

RGiL±,Ly)(xL2 , (1) 
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hrms{L±,Ly) = Llh((ALy)/{AL^y) , 



where the universal exponents z 
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< 1, A is an ultraviolet cutoff, /^(m) constant for 
M — > and //i(w) oc u^"^, for m — > oo. 

For general spatial dimension d, Flory theory treat- 
ment of the self-avoidance (SA) predicts vi = sug- 
gesting that the tubule phase should be stable down to 
the lower critical dimension die = 2, where i/ = 1, and 
therefore should exist in 3d, predicting = which 
implies that 2; = |, and C = However, the analo- 
gous Flory result for the crumpled phase v = has 
more recently been called into question. Numerical sim- 
ulations [6] find no crumpled phase below d = 4. An 
uncontrolled Gaussian approximation [7,8] supports this 
finding, predicting v = 2^ which suggests that du = 4 for 
the crumpled phase. Both this, and the numerical values 
of V for d > 4, agree well with the simulations. 

The same approximation for the tubule phase [4] gives 
v = , which, unfortunately, suggests that die = 4 
for the tubule phase as well. However, despite its success 
for the crumpled phase, this Gaussian approximation is 
known to be far from trustworthy. For example, it pre- 
dicts v = 2 for linear polymers, which not only is less 
accurate than the Flory result v = but also incor- 
rectly predicts die = 2 for polymers, when, in fact, it is 
well known that die = 1 in that case. 

Whether the Gaussian approximation is any more re- 
liable for our tubule phase remains an open question. 
One could argue that a slice perpendicular to the y-axis 
through our tubule looks like a SA random walk in two di- 
mensions, for which the Flory result oi v = 3/4 is known 
to be exact, while the Gaussian approximation = 1 is 
clearly wrong. Whether or not this analogy holds, it is 
clear that SA, though a relevant perturbation, is much 
less important for tubules than for the crumpled phase, 
since points on the membrane widely seperated in the y- 
direction never bump into each other in the tubule phase, 
while they do in the crumpled phase. So it seems quite 
plausible, the Gaussian approximation notwithstanding, 
that the tubule phase is stable in d = 3. Furthermore, 
the suppression of the crumpled phase by SA, makes the 
possibility of the new tubule phase even more interesting 
and important. Whether the tubule phase does survive 
in d = 3 can only be determined by simulations and ex- 
periments on anisotropic membranes, both of which we 
hope our work stimulates. 

In the following discussions, numerical estimates for 
the values of the exponents will be obtained from the 
Flory estimate z^j = | in 3d; these numbers should be 
taken with a grain of salt, due to the uncertainties just 
discussed, about the validity of the Flory theory. 

Eq.2 implies that hrms{L) oc L^""' « L^/^, for a 



square membrane with Lj 



L, using 



ALy/ {AL±J^ — > 00. In the "linear polymer" limit 

T T , h nr ,.z(3/2-C) 
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Ly^^ /L^/'^ . Defining Lp to be the value of Ly at which 
hrms = Ly,we obtain the orientational persistence length 
Lp oc L'^_^'~^^ « L^^^. For any roughly square membrane, 
Ly is much less than Lp, hence the tubule phase is stable 
against thermal fluctuations as L — )■ 00. 

Like the flat phase [1,2,9], the tubule phase exhibits 
anomalous elasticity; however, as discussed above, SA is 
a strongly relevant perturbation in this new phase. The 
tubule, swelled by the SA interaction, acts for Ly >> 
Ly = A~^(ALj_)^ Hke a polymer with bending rigidity 

Kp(L_l) oc L_l-Rg « -f'x"^''' ~ ^^r- 

Our model for anisotropic membranes is a generaliza- 
tion of the isotropic model [11]. We characterize the con- 
figuration of the membrane by the position f(x) in the 
d-dimensional embedding space of the point in the mem- 
brane labelled by a D-dimensional internal co-ordinate 
X. In the physical case, D = 2 and d = 3. The Lan- 
dau free energy F is an expansion in the local tangent 
vectors daf{Tc), keeping only the leading terms consistent 
with global translation and rotation invariance: 



m^)] = \j d^'-'x^ 
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+ ^ {d^r- d^fY + v^y {d^fY {dyr-f 



Ifd'^xj d^x' <5('^)(r(x) -r(x')) 
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where the k's, t's, u's, u's are elastic constants, and h is 
the SA interaction strength. The first three (/t) terms in 
F represent the anisotropic bending energy of the mem- 
brane. The clastic constants t±_ and ty are > at high 
temperatures and < at low temperatures. When both 
are positive, the membrane crumples. When either is 
negative, the membrane extends in the associated direc- 
tion, and the u and v quartic terms are then needed 
for stability. Eq.3 reduces to the model for isotropic 



membranes [11] when t±_ = ty, k^±_ 



K±y 



0, 



4(w -I- u), u±± = u±y = Au, and = v^y 



Av. 



In mean-field theory, we seek a configuration r (x) that 
minimizes the free energy Eq.3. Since the curvature en- 
ergies /t_L {d^r) and Ky (dyf) vanish when r(x) is lin- 
ear in X, we seek these minima by inserting the ansatz 
r^(x) = {(±x±,(yy,0,0, . . . ,0) into Eq.3. For now ne- 
glecting the SA interaction. 



-L^ ^Ly 



tyCl+t^{D-l)Cl 



+ ^<±iD - ifCl + ^UyyC'y + V ^y{D - 1 ) ^ 



• (4) 
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where u'. 



v±± + u±±/{D — 1). Minimizing F over 



and (y yields two possible phase diagram topologies. For 
u'j_j_Uyy > v']_y, we obtain Fig.l. Both and (y vanish 
for t±, ty > 0. This is the crumpled phase: the entire 
membrane, in mean-field theory, collapses into the origin, 
= = i.e., r(x) = for all x. In our new y-tubule 
phase, characterized by = and Q = \/\ty\/uijy > 0: 
the membrane is extended in the y-direction but crum- 
pled in the _L directions. The _L-tubule phase is the anal- 
ogous phase with the y and _L directions reversed. The 
tubule-flat boundary slopes are Uyy/v±y and v±y/u'j_j_ 
respectively. In the flat phase, both and (y ^ 0. 
For u'j_j_Uyy < v]_y, the flat phase disappears, and is re- 
placed by a direct flrst-order transition from _L-tubule to 
i/-tubule along the locus ty = {v±y/u'j_j^t±. 

The flat and the crumpled phases of anisotropic mem- 
branes in Fig.l are in the same universality class [10] 
as those of isotropic membranes [1]. In the crumpled 
phase, t±, ty > 0, and all other local terms in Eq.3 are 
irrelevant at long wavelengths. A change of variables 
xj_ = x'^ytj_/ty makes the remaining energy isotropic. 

We now consider the effects of fluctuations, ignor- 
ing SA (i.e., the "phantom" membrane). Consider 
the y-tubule phase. To treat fluctuations, we perturb 
around the mean-fleld solution by writing ^(x) = [CyV + 
m(x),/i(x)), where h{x) is a d — 1-component vector 
orthogonal to y. Inserting the above expression for f 
into Eq.3 and keeping only relevant terms gives Ftot = 
Fmft + F^i, where F^ft = \Ll^^Ly[tyCl + jUyyC^], 



x±dy 



+ gy{dyu + -{dyhff 



(5) 



K = Ky, t = t^ + V^yCl, Qy = %yCy/2, aUd g_L = 
t + U±yCfj. 

Note that the ratios of the coefficients of the quadratic 
(dyu)'^ and the anharmonic dyu{dyh)'^ and {dyhY terms 
in Fei must be exactly 4:4: 1, since they must appear 
together as a result of expanding the rotationally invari- 
ant combination [dyU + ^idyhY)"^ . This ratio aflows us to 
calculate exactly the long-wavelength anomalous elastic- 
ity of "phantom" tubules, as we wfll show in a moment. 

The propagators can be read off from Eq.5, giv- 
ing (/l,:(q)ftj(-q)) = kBT5iPh{(\), (u(q)u(-q)) = 



fcBTG„(q), where G^i(q) 



g^_q\ + gyQy, and Sjj is a Kronecker delta when both 
indices i and j ^ y, and is zero if either i or j = y. The 
rms fluctuations in the harmonic approximation are: 



/ 

Jl± 



d q±dqy 



{2Tr)D tql + Kql 



4 °^ ^1. 



(6) 



clearly revealing that the upper critical dimension is 
Due = 5/2. Below D^c, we expect anomalous elasticity. 



However, this anomaly is not manifested in the fluctua- 
tions of h alone. We can see this by integrating out the 
phonons u exactly, the only remaining anharmonicity in 
the effective elastic free energy for h alone is, 

Fanh[h] = \ I (K(ki) ■ K(k2)) fK(k3) • h{ln)] X 

X (ki ■ ka) (kg ■ k4) y^q) , (7) 

where q = ki -I- k2 and ki + k2 -f ka -I- k4 = 0. The 
effective vertex is V'ft(q) = gyg±q]_l {gyql + g±q\), and is 
irrelevant for D > 3/2, as can be seen by simple power 
counting. Thus, in D = 2, the elastic constants t, g± and 
Ky are finite and non-zero as ^ 0. 

However, gy is driven to zero as qy 0. In a self- 
consistent one-loop perturbative calculation, similar to 
that successfully used to compute the anomalous elastic- 
ity in the fiat phase [9], we find 



fly(q) = gy 



(8) 



/ 



{kBTrgl{ci)pl{Py - qy? d^-'p^dpy /{2Tr)^ 
'tpl + k(p>4) {t\ps. - q±P + Ki\p- q\)ipy - qyY) ' 



where g° is the "bare" value of gy. The above argument 
shows that K,{p) can be replaced by a constant, since the 
h elasticity is not anomalous. This self-consistent equa- 
tion can be solved by the ansatz, gy(q) = Qy" ^g(qy/q±)■ 
with z = ^, rju = 5 — 2D, which we have verified works 
to all orders in perturbation theory. 

We now compute the phantom tubule diameter Rq and 
transverse wandering roughness hrms, defined by Rq = 
{\h(L^,y)-hiO^,y)['), hl^, = (|/I(x±,L,)-Mxx,0)h- 
Because hrms and Rq receive large contributions from 
q_L = and f/y = zero modes, respectively, Rg and 
hrms , surprisingly, scale in different ways with the mem- 
brane dimensions L± and Ly. Taking into account the 
zero modes, we calculate hrms and Rq by equiparti- 
tion, and find the forms Eqs.1,2, with z = 1/2, i/t = 
(5 — 2D) I A, and C, = 2i/t. For a nearly square mem- 
brane Lj_ ^ Ly ^ L ^ 00, for which hLy >> (ALj^)^, 



we obtain hrms oc Lyi'^f ^^^^ oc L^--^/^ Thus for a 
D = 2 phantom tubule, hrms c< L. UnHke the fiat phase, 
no log(L/a) correction arises, so the {D = 2) phantom 
tubule is just marginally stable, but with wild transverse 
undulations. These arc greatly suppressed by SA, to the 
effects of which we now turn. 

We begin by estimating the radius of gyration using 
Flory theory. Specializing henceforth to the physical 
case £) = 2 for simpHcity, we estimate the SA energy 
Eg A as EsA ~ bVp"^ , where the volume V in the embed- 
ding space occupied by the tubule is F « Rff^Ly, and 
the density p in the embedding space of the tubule is 
p « M/V oc L±Ly/V, where M is the mass of the mem- 
brane. Putting these together gives Esa oc LyL']_/R'^^. 
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Inserting our earlier, phantom membrane result Rg oc 
(for D = 2) and taking L±_ oc Ly, as required 

by anisotropic scaling, we find that Esa oc 
which goes to infinity as Ly ^ oo for d < 11. Thus, 
SA is strongly relevant, and changes the long wavelength 
behavior of the membrane, for d < 11. 

We can calculate Rq for d < 11 by combining the 
above estimate of Esa with a similar scaling estimate of 
the elastic energy, yielding: 

= [tyCl + UyyC'y + t{^)']L^Ly + . (9) 

Minimizing this over Rq, we obtain Rg{L±) oc L'^, with 
Ut = For the physical case d = 3, this gives Ut = 3/4. 
Since a slice through the tubule traces out a crumpled 
polymer embedded in 2d, we conjecture that Ut = 3/4 
is an exact result for the tubule thickness, since it is for 
2d polymers. For a square membrane, Ly ~ L±, it is 
straightforward to argue that the Qy = zero modes do 
not contribute to Rq, and L± is the dominant infra-red 
cutoff. Hence Eq.l gives the correct radius of gyration. 

This highly nontrivial ground state for the SA tubule 
is not modified by thermal fiuctuations. That is, even at 
T > 0, the variation of Rq with L± is completely domi- 
nated by the SA energy, i.e. determined by a T = fixed 
point. This can be seen by evaluating the elastic (or SA) 
energy with Rg{L±) given by this SA ground state. For 
an L X L, membrane, this energy Egi gc , 9 = = | 
(for d = 3), is much larger than fc^T as L — > oo. 

Using a generalization of Landau's derivation of 
shell theory we now calculate anomalous elasticity 
in the presence of SA. Bending the tubule with ra- 
dius of curvature Rc 3> Rg, induces a strain e ~ 
Rg/Rc, which costs an additional elastic energy den- 
sity gy{Ly,L±)e'^ =gy{Ly,L±){RG{Ly)/Rc)'^. Interpret- 
ing this additional energy as an effective bending energy 
density Ky{L±, Ly) / R^ leads to the effective bend mod- 
ulus Ky{L±, Ly), 

Ky{Lj_,Ly) ~ gy{LA_,Ly)RG{Lj_,Ly)'^ , (lO) 

Inserting Ky{L±,Ly) = L'^f^{ALy/{AL±)^), 
gy{L^,Ly) = L-^-fg{KLyl{kLj_)^) and RG{L±,Ly) = 
L'2fR{ALy/{AL±)^) into above expression, we obtain a 
relation between the exponents 2f( = z{7] + 7]u), which is 
satisfied by our earlier results for the phantom tubule. 

The physical SA tubule at T = is absolutely straight 
(i.e. = 1 exactly). This impHes that the tubule stretch- 
ing elastic constant Qy = Qy, its bare value, since there are 
neither fluctuations nor SA effects to renormalize it (in 
contrast to k). Therefore = at T = 0; Eq. 10 then 
implies that Kq for the T = SA tubule is already anoma- 
lous, and given by /to(q) = qy'^"'^'' f^{qyA^ /{Aql)) . 

The effective free energy describing thermal fiuctua- 
tions about this nontrivial, SA ground state is 




Balancing the Kqy term in this expression with the tq"]^ 
term gives the anisotropy of scaling exponent z, defined 
by qy oc q^. We thereby obtain z = 2-lt/z ' '^hich 
gives z = Now calculating the fiuctuation cor- 

rections to gy from Eq.8 using the wavevector dependent 
K found above, we find that the integral in Eq.8 con- 
verges in the infra-red and Qy is finite as q pro- 
vided 2 > |(1 — z/j). If is anywhere near its Flory 
value = 3/4 in D = 2, d = 3, this condition is clearly 
satisfied. Thus, unlike the phantom, the SA tubule has 
riu = 0, even at T > 0. Using this fact in Eq.lO and 
taking the Flory expression for ut in d = 3, we obtain. 

Using Eq.ll to compute the transverse tubule undu- 
lations hrms, we find hrms = L'^fh{^Ly/{AL±)^), where 
C = I - and Mx) « a;(i+2'")/2^ and h{x) « con- 

stant for a; —> 00 and a; — > 0, respectively. For a roughly 
square L x L membrane, we are always in the a; — > oo 
limit of the crossover function fh{x), which implies that 
the rms transverse undulations of the tubule are given 
by hrms ~ L3/2-(i+2^0/2 ^ ^1/4, Similar arguments 
applied to Rg yield Eq.l. 

We [4] are currently investigating the effects of SA on 
the tubule phase in d = d^f — e, following the work of 
Ref. [12] for the crumpled phase. The scaling theory of 
the tubule to fiat and crumpled to tubule transitions, 
along with a Flory theory for these transitions, will be 
presented in a future publication [4]. 
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